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Abstract 

Background: Laminopathies are diseases characterized by defects in nuclear envelope structure. A well-known 
exannple is Emery-Dreifuss muscular dystrophy, which is caused by mutations in the human lamin A/C and emerin 
genes. While most nuclear envelope proteins are ubiquitously expressed, laminopathies often affect only a subset 
of tissues. The molecular mechanisms underlying these tissue-specific manifestations remain elusive. We 
hypothesize that different functional subclasses of genes might be differentially affected by defects in specific 
nuclear envelope components. 

Results: Here we determine genome-wide DNA association profiles of two nuclear envelope components, lamin/ 
LMN-l and emerin/EMR-l in adult Caenorhabditis elegans. Although both proteins bind to transcriptionally inactive 
regions of the genome, EMR-l is enriched at genes involved in muscle and neuronal function. Deletion of either 
EMR-l or LE)\/l-2, another integral envelope protein, causes local changes in nuclear architecture as evidenced by 
altered association between DNA and LMN-l. Transcriptome analyses reveal that EMR-l and LEM-2 are associated 
with gene repression, particularly of genes implicated in muscle and nervous system function. We demonstrate that 
emr-1, but not lem-2, mutants are sensitive to the cholinesterase inhibitor aldicarb, indicating altered activity at 
neuromuscular junctions. 

Conclusions: We identify a class of elements that bind EMR-l but do not associate with LMN-l, and these are 
enriched for muscle and neuronal genes. Our data support a redundant function of EMR-l and LEM-2 in chromatin 
anchoring to the nuclear envelope and gene repression. We demonstrate a specific role of EMR-l in neuromuscular 
junction activity that may contribute to Emery-Dreifuss muscular dystrophy in humans. 



Background 

The nuclear envelope (NE) is an essential structure con- 
served in all eukaryotes that separates the nucleus from 
the cytoplasm [1], It is characterized by the presence of 
two concentric lipid membranes, the outer nuclear mem- 
brane (ONM), which is continuous with the endoplasmic 
reticulum, and the inner nuclear membrane (INM). The 
two membranes are fused where large nuclear pore com- 
plexes create transport channels to control trafficking of 
molecules between the nucleus and the cytoplasm. Numer- 
ous membrane proteins are integrated in the membrane 
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bilayer. In the ONM, membrane proteins interact with mi- 
crotubules, actin filaments and other cytoplasmic compo- 
nents. In the INM of higher eukaryotes, membrane 
proteins interact with the nuclear lamina, a protein mesh- 
work of intermediate filaments that serves as a scaffold to 
maintain the morphology of the nucleus. INM and nuclear 
lamina proteins also interact with chromatin [2]. 

The NE is an important player in a number of cellular 
functions, including maintenance of structural integrity 
of the nucleus, chromatin organization, transcriptional 
regulation, DNA replication, cell-cycle control, differen- 
tiation, nuclear migration, and apoptosis [1,3-5]. Reflect- 
ing these many functions, mutations in NE components 
have been associated with numerous human diseases. In 
particular, a predominant group of diseases linked to 
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alterations in the NE is caused by mutations in genes 
encoding the nuclear lamins {LMNA, LMNBl and LMNB2 
in humans) and is referred to as laminopathies. Lami- 
nopathies include diverse pathologies, such as muscular 
dystrophies, cardiomyopathies, bone disorders, neuro- 
logical diseases, lipodystrophies and premature aging 
[5,6]. Laminopathy mutations have also been identified 
in several genes encoding ONM and INM proteins. 
Emery-Dreifuss muscular dystrophy (EDMD) was initially 
linked to mutations in the EMD gene, which encodes the 
INM protein emerin [7], but can also be caused by muta- 
tions in LMNA and NE genes SYNEl, SYNE2 and 
TMEM43 [6]. EDMD develops during childhood and 
manifests as degenerative muscle atrophy with cardiac de- 
fects [6]. Although clinical investigations focus on muscle 
abnormalities in EDMD patients, mice lacking emerin 
were found to display minimal muscle dysfunction but di- 
minished motor coordination [8], suggesting that emerin 
could have additional functions apart from regulation of 
muscle formation and integrity. Interestingly, differences 
between humans and mice in relative expression levels of 
emerin and lamina-associated polypeptide 1 (LAPl) were 
recently proposed to explain the absence of muscle dys- 
trophy in emerin mutant mice [9]. Since most NE genes 
are expressed ubiquitously, a main challenge is to under- 
stand how laminopathy mutations result in tissue-specific 
phenotypes [5,6]. One proposed hypothesis argues that an 
alteration of the integrity of the NE - for example, by af- 
fecting the nuclear lamina or the connection of the NE 
with the cytoskeleton - could lead to a structural weakness 
and a decrease in the ability of the nucleus to resist mech- 
anical stress [10-12]. This could be of special relevance in 
muscle and skeletal disorders. However, symptoms affect- 
ing others tissues, such as the nervous system, cannot be 
explained simply by mechanical problems. Instead, several 
studies indicate that the NE contributes to chromatin 
organization by direct interaction with chromatin factors. 
For instance, the INM LEM-domain-containing proteins 
(LAP2, emerin and MANl) interact with the chromatin 
binding protein BAP (also known as BANFl), histone dea- 
cetylases as well as with the nuclear lamins [2,13,14]. The 
nuclear lamins also interact with chromatin remodelers 
involved in heterochromatin formation, such as hetero- 
chromatin protein HPl or the NURD complex [15]. 

Analyses of different organisms and cell types have 
shown that the NE is mainly associated with hetero- 
chromatin and silenced DNA [16-19]. Consistent with 
this, alterations in the distribution of heterochromatin 
marks, such as histone H3 lysine 9 methylation (H3K9me) 
and H3K27me, have been observed in murine and hu- 
man myoblast cells harboring mutations in the LMNA 
gene [20,21]. Mutations in NE proteins can also lead to 
a reduction of transcription of muscle differentiation 
genes [22]. Although this reduction could at least 



partially be due to the role of NE components in het- 
erochromatin formation, it has also been proposed that 
INM proteins can sequester transcription factors to the 
nuclear periphery and impede their binding to target 
genes. For example, in humans and mice, emerin phys- 
ically interacts with lmo7 and p-catenin, two transcrip- 
tion factors involved in muscle differentiation [23,24], 
whereas in humans LEMD3/MAN1 tethers Smads to 
the NE, thereby affecting connective tissue differenti- 
ation [25-27]. 

While the importance of the NE as a regulator of nu- 
clear architecture and gene expression is becoming in- 
creasingly evident, the molecular mechanisms underlying 
this remain elusive. In particular, the contributions of indi- 
vidual NE proteins to tissue-specific functions are not well 
understood. Characterization of the chromatin domains 
that interact with the NE is required to decipher how the 
NE contributes to nuclear organization. Until now, most 
experiments have focused on cultured cells, whereas few 
studies have been performed on cells within intact organ- 
isms. We therefore decided to analyze the DNA associated 
with two components of the NE, lamin/LMN-1 and 
emerin/EMR-1, in whole adult Caenorhabditis elegans. 
Due to its amenability to genetic manipulations, C. elegans 
is particularly suitable to genomic analyses across several 
genotypes and developmental stages [28]. Our results 
show that both LMN-1 and EMR-1 are associated with 
lowly expressed genes. As expected, similar DNA profiles 
were observed for the two proteins, but we also identified 
elements bound by only one of the two. 'EMR-1 only ele- 
ments were enriched for muscle and neuronal genes, 
which became accessible to LMN-1 association when 
emr-1 was deleted. Furthermore, we observed that EMR-1 
acts redundantly with another LEM domain containing 
protein, LEM-2, to repress transcription, consistent with 
their functional redundancy during mitosis [29], develop- 
ment and myogenesis [30], and signaling [31]. Finally, we 
demonstrate that EMR-1, but not LEM-2, is required for 
proper neuromuscular junction activity. Together, this 
study demonstrates the importance of EMR-1 in the con- 
trol of chromatin organization and gene expression of 
muscle and neuronal genes, thereby providing clues as to 
how defects in INM protein function can have tissue- 
specific consequences. 

Results 

Identification of chromatin anchored to LMN-1 and EMR-1 
by DamID 

To investigate the specific role that different components 
of the NE play in the control of chromatin organization 
and gene expression, we generated genome-wide inter- 
action maps of the sole C. elegans nuclear lamina protein, 
lamin/LMN-1, and the inner nuclear membrane protein 
emerin/EMR-1, in adult C. elegans nematodes using 
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DamID [32]. This technique is based on the in vivo ex- 
pression of chimeras of the Escherichia coli DNA adenine 
methyltransferase (Dam) and a chromatin-interacting pro- 
tein. Upon interaction of the Dam-fused protein with 
chromatin, adenines in the vicinity are methylated. These 
DNA regions are subsequently identified by microarray 
analysis or high-throughput sequencing. DamID has been 
successfully used in many organisms and has been dem- 
onstrated to reliably identify NE-associated sequences 
[17,19]. To minimize experimental variation, we created 
C elegans strains containing single copy insertions of the 
chimeric transgenes in chromosome II (Figure SI A in 
Additional file 1; Materials and methods). We used the 
promoter of the heat shock-inducible gene hsp- 16.41 , and 
due to the requirement for only trace expression levels of 
the chimeric protein, performed the DamID experiments 
without induction [19]. Without heat shock, transgene ex- 
pression is undetectable by immunofluorescence or west- 
ern blot analyses: however, after heat shock, the nuclear 
peripheral localization of chimeric proteins was readily 
observed (Figure SIB, C in Additional file 1). We also gen- 
erated a control strain expressing GFP::Dam (which freely 
diffuses throughout the nucleus; Figure SID in Additional 
file 1) to normalize for local chromatin accessibility [19]. 
Under normal growth conditions, nematodes containing 
chimeric proteins were fertile and indistinguishable from 
wild-type (Bristol N2) animals. We conclude that the 
Dam-fusion proteins mimic the localization of their en- 
dogenous counterparts and their presence did not pro- 
duce any apparent deleterious effects. 

To investigate nuclear organization in differentiated 
cells we performed DamID on young adults. Thus, the 
DamID signals we describe in the sections below repre- 
sent the mixture of cell types of which the nematode 
is composed. We isolated adenine-methylated DNA 
from three biological replicates for each strain (Dam:: 
LMN-1, Dam::EMR-l and GFP::Dam) and analyzed it on 
C. elegans whole-genome high-density tiling microarrays. 
We normalized the DamID signals using MA2C [33] 
(Materials and methods) and then all MA2C data were 
quantile normalized to facilitate comparison between 
strains. The replicas were highly correlated (Pearson cor- 
relation between replicates ranging from 0.88 to 0.98; 
Figure S2A in Additional file 1) and genomic profiles of 
LMN-1 and EMR-1 association were constructed by 
averaging the three replicates. 

Maps of LMN-1 and EMR-1 reveal global similarities and 
association with silent DNA 

The genome-wide profiles of Dam::LMN-l and Dam:: 
EMR-1 signals were very similar to each other (Pearson 
correlation of 0.88; Figure lA). Both LMN-1 and EMR-1 
associate preferentially with autosomal arms within 4 to 
7 Mb of the telomeres. In contrast, the central part of 



the autosomes was largely devoid of LMN-1 and EMR-1 
association (Figure IB; Figure S3A in Additional file 1). 
Interestingly, the X chromosome showed a different pat- 
tern since interaction was mainly observed at the left 
end of the chromosome whereas only weak association 
was seen at the right arm. These patterns were similar to 
those of INM proteins LEM-2 [18] and LMN-1 [34] re- 
ported in embryos. Thus, the organization of autosomes 
within the nucleus is fundamentally different from that 
of the X chromosome in C. elegans. 

To precisely define LMN-1 -associated domains (LADs) 
and EMR-1 -associated domains (EADs), we converted 
positive and negative DamID scores to +1 and -1, respect- 
ively, and used a window-based method to identify DNA 
regions continuously associated with LMN-1 or EMR-1 
[18] (Materials and methods). Using a false discovery 
rate <5%, we identified 1,021 LADs and 1,106 EADs, each 
of which covers 37.5% and 37.8% of the genome, respect- 
ively (Figure IC; Table SI in Additional file 2). The size of 
LADs and EADs ranged from 8 to 300 kb with a median 
size of 22.7 kb and 23.9 kb for LMN-1 and EMR-1, re- 
spectively (Figure ID). A large proportion of LADs and 
EADs was overlapped: 89% of DNA in LADs was also as- 
sociated with EMR-1. LADs and EADs at the chromo- 
some arms were separated by non-associated regions 
('gaps'), which have a median size of 13.5 kb and 19.2 kb 
for LMN-1 and EMR-1, respectively (Figure ID). It is im- 
portant to note that genomic features, such as open read- 
ing frames, promoters and intergenic regions are found in 
both domains and gaps (Figure IC). 

Like LEM-2-associated regions [18], the fraction of 
DNA regions occupied by LMN-1 or EMR-1 in auto- 
somes was positively correlated with the size of the auto- 
somes (Figure IE). In addition, we observed an equal or 
higher proportion of each autosome bound by LMN-1 
than EMR-1. In contrast, the X chromosome showed a 
lower association with LMN-1 compared to EMR-1, sug- 
gesting that EMR-1 might play a distinct role in anchor- 
ing the X chromosome to the NE (Figure IE; see below). 

Previous studies in C. elegans embryos and other or- 
ganisms have shown that the chromatin regions associ- 
ated with the NE are marked predominantly by histone 
modifications characteristic of heterochromatin and are 
transcriptionally silent [17-19]. To investigate if this dis- 
tribution is maintained in adult nematodes, we analyzed 
the levels of H3K4me3 and H3K36me3, both associated 
with transcriptional activity [35], in LADs and EADs. 
Sliding window analyses across EMR-1 domain boundar- 
ies showed that both LADs and EADs have lower 
H3K4me3 and H3K36me3 levels, and instead, gaps have 
higher levels of these two marks, suggesting that tran- 
scriptionally active chromatin is mainly localized in the 
nuclear interior in C. elegans adults (Figure IF; Figure 
S3B in Additional file 1), as it is in embryos [18]. 
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Figure 1 LMN-I and EMR-I have overlapping DNA association profiles devoid of active chromatin marks. (A) Genome-wide correlation 
plot of LMN-1 and EMR-1 occupancy. The correlation coefficient (r) is at the lower right corner. (B) Regions enriched for Dam::LMN-l and Dam:: 
EMR-1 in chromosomes V and X. For each track, the average MA2C score probe signal of three independent biological replicates is plotted. Other 
chromosomes are shown in Figure S3 in Additional file 1. (C) Representative region enriched for LMN-1 (MA2C score; blue) and EMR-1 (MA2C 
score; green) association. Positions of annotated genes are indicated below. LMN-1 and EMR-1 -associated domains (LADs and EADs, respectively; 
black bars) were defined by a sliding window algorithm. (D) Sizes of LADs, EADs and regions between domains (Gaps). Bottom and top of boxes 
indicate the 25^^ and 75^^ percentiles, respectively, and lines in the boxes indicate medians. Whiskers indicate the lowest or the highest data point 
within 1.5x interquartile range from the box. Probability value from Wilcoxon rank sum test is indicated: **P< 10"^. (E) Relationship between 
chromosome size and LMN-1 or EMR-1 relative occupancy (sum of DamID domains divided by chromosome size). Linear regression lines for 
LMN-1 (blue; r = 0.82) and EMR-1 (green; r = 0.93) occupancy on autosomes are indicated. (F) Average H3K4me3 and H3K36me3 profiles at EAD 
boundaries. Chromatin immunoprecipitation (ChlP)-chip probe signal scores for indicated histone modifications are plotted. For comparison, 
averages of LMN-1 and EMR-1 DamID MA2C score are also shown. 
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However, low H3K4me3 and H3K36me3 levels were 
not exclusively observed in LADs and EADs but also 
inside gaps, suggesting that not all chromatin away 
from the NE is transcriptionally active (Figure S3B in 
Additional file 1). 

Chromosome architecture is maintained during 
development 

Given the apparent similarity of chromosome organization 
in C. elegans embryos and adults, we analyzed in more de- 
tail the pattern of chromosomal associations with the NE 
in these two developmental stages. First, we compared the 
embryonic chromatin association profile of LMN-1 ob- 
tained by DamID [34] with the one of LEM-2 generated 
by chromatin immunoprecipitation (ChIP)-chip/seq [18]. 
The highly similar profiles cross-validate the two tech- 
niques (Figure 2A). Second, we also noted that the global 
profiles observed by LMN-1 and EMR-1 DamID during 
adulthood were very similar (Figure 2B), suggesting that 
the two DamID fusion proteins reliably detect chromatin 
at the NE. Finally, we compared LMN-1 DamID profiles 
in embryos and adults. In general, the boundaries between 
LMN-1 -associated arms and non-associated centers coin- 
cided well in the two developmental stages. In addition, 
the dynamic range of LMN-1 association at chromosome 
arms was very similar between the two developmental 
stages (Figure 2C, D). However, in adulthood we observed 
slightly more NE association in the arms of chromosome 
IV, and during early embryogenesis higher association in 
the arms of chromosomes I, II and III (Figure 2C, D). 

Interestingly, we found a clear difference in the behav- 
ior of the central part of the chromosomes. In early em- 
bryos, LMN-1 MA2C scores for the central part of the 
chromosomes were close to zero, while in adult nema- 
todes LMN-1 association with chromosome centers was 
significantly lower (Figure 2C, D). Taken together, these 
data suggest that during adulthood segregation of chro- 
matin into NE-associated and NE-excluded domains is 
more pronounced, especially at the central region of the 
chromosomes, which seems to be excluded from the nu- 
clear periphery in most cells. 

Local differences between LMN-1 and EMR-1 DNA 
association reveal a specific role of EMR-1 in the 
regulation of muscle and neuronal genes 

Despite the global similarities between LADs and EADs, 
there were also differences in the DamID profiles during 
adulthood. Specifically, the DamID signal of LMN-1 was 
higher than that of EMR-1 on the left arm of chromo- 
somes I and IV and on both arms of chromosome III 
within regions covering 400 kb (Figure 3A). To perform 
a more detailed analysis of the differences in the chro- 
matin association of LMN-1 and EMR-1, we subtracted 
the MA2C values obtained with the two Dam fusion 



proteins and plotted the resulting values (Figure S4A in 
Additional file 1). This revealed that areas of preferential 
association with LMN-1 were strongly enriched at auto- 
some arms and at the left end of the X chromosome. 
Most likely, anchoring of these regions to the NE is me- 
diated primarily by direct interaction with LMN-1 or by 
association with NE proteins other than EMR-1. 

Although LMN-1 is ubiquitously expressed and phys- 
ically interacts with EMR-1 [36], we also found elements 
associated exclusively with EMR-1 (Figure 3B, C). Con- 
trary to LMN-1 -enriched regions described above, these 
elements were more abundant close to the central part 
of the autosomes and on the X chromosome (Figure 3B; 
Figure S4 in Additional file 1). A systematic search 
identified 340 'EMR-1 only elements (Materials and 
methods) ranging from 339 bp to 3,328 bp and with a 
median size of 763 bp (Table SI in Additional file 2). 
These elements were included in EADs and surrounded 
by LADs (Figure 3C; Figure S4B in Additional file 1). 
'EMR-1 only' elements occurred significantly more often 
in promoters than expected, and occurred less often 
than expected in exons and introns (Figure 3C, D). 
Interestingly, Gene Ontology (GO) analyses of genes 
placed in or immediately downstream of EADs identified 
29 GO categories, most of which were related to muscle 
and neuronal physiology, including signaling transduc- 
tion at synapses (G-protein coupled receptor protein 
signaling pathway, P = 2.41e'^; synaptic transmission, 
P = 0.01; acetylcholine catabolic process, P = 0.04), be- 
havior (behavior, P = 0.04), reproduction (reproductive 
process, P = 0.01) or locomotion (regulation of locomo- 
tion, P = 0.006) (Figure 3E; Table S5 in Additional file 3). 
Employing the same criteria used to defined 'EMR-1 
only elements, we identified 97 'LMN-1 only' elements 
ranging from 406 bp to 2,097 bp and with a median size 
of 661 bp (Table SI in Additional file 2). These 'LMN-1 
only' elements were located almost exclusively in the large 
regions with high LMN-1 /EMR-1 DamID ratios previ- 
ously described (Figure 3B; Figure S4C in Additional 
file 1) and they were not enriched for genome features or 
GO term categories. 

The identification of 'EMR-1 only' elements indicates 
that EMR-1 could have a predominant function in certain 
tissues. Indeed, we have observed that EMR-1 expression 
varies significantly between tissues (Figure S5 in Additional 
file 1). To further investigate whether LMN-1 and EMR-1 
could play tissue-specific roles in chromatin organization, 
we analyzed LMN-1 and EMR-1 DamID scores of 6,954 
genes that are expressed in a tissue-specific manner, as 
previously defined by others (Table S2 in Additional file 2; 
see Materials and methods for references). Of 6,954 genes, 
4,909 (70.6%) were not anchored to the NE as reflected by 
the low MA2C values derived from DamID experiments 
performed in the whole animal (Figure S6A in Additional 



Gonzalez-Aguilera et al. Genome Biology 2014, 15:R21 
http://genonnebiology.conn/201 4/1 5/2/R21 



Page 7 of 1 9 



oT 1.1, 

£^ 0 

CD O ^ 
LU < 

e.-i.6 



LMN-1 LEM-2 



B 



■o 1.4 
1..1.5 



Chrl Chrll Chrlll ChrlV 
LMN-1 EMR-1 



ChrV 



ChrX 




Chrl Chrll Chrlll ChrlV 
Adult Embryo 



ChrV 



ChrX 



Chrl Chrll Chrlll ChrlV 



ChrV 



ChrX 



D 



o 





■^HiT ■■r' -i"^ 







EA AA EC AC EA AA EC AC EA AA EC AC EA AA EC AC EA AA EC AC EL AL ER AR 



Chrl 



Chrll 



Chrlll 



ChrlV 



ChrV 



ChrX 



Figure 2 Conservation of nuclear organization during development. (A, B) Global association with chromatin is highly similar among NE 
proteins, both in embryos (A) for LMN-1 (black) and LEM-2 (orange) and during adulthood (B) for LMN-1 (blue) and EMR-1 (green). Embryonic 
LMN-1 DamID and LEM-2 ChIP MA2C scores were obtained from [34] and [18], respectively. [C, D) Representation of LMN-1 association as 
chromosomal profiles (C; embryos in black, adults in blue) and boxplot (D; as in Figure 1). (D) LMN-1 MA2C scores in embryos (E-) and adults (A-) 
averaged over chromosome arms (-A) and centers (-C), or, for the X chromosome, over the left (-L) and right (-R) portions. Probability values 
from Wilcoxon rank sum tests are indicated: *P< 0.05, ***p< 10"". DamID signals at chromosome ends do not differ dramatically between the 
two developmental stages, whereas chromosome centers are devoid of LMN-1 interaction, particularly in adults. Probe scores were averaged in 
100-kb windows for all panels. 



file 1). Importantly, in the 2,045 genes with a positive 
MA2C score (that is, anchored to the NE in most cells of 
the nematode), we observed significantly more association 
of intestinal, germline, and ubiquitous genes with LMN-1 
than with EMR-1 (Figure 3F). In contrast, muscle and 
neuronal genes were associated with EMR-1 to a higher 
degree than they were with LMN-1 (Figure 3F; Figure S6B 
in Additional file 1), which is consistent with the GO ana- 
lysis of 'EMR-1 only elements. Our data thus suggest that 
EMR-1 is specifically involved in gene regulation in muscle 
and neurons. 



Lack of EMR-1 or LEM-2 produces local changes In nuclear 
organization 

The unexpected finding of 'EMR-1 only elements led us 
to speculate how they can be established. We imagine 
three scenarios for their formation: (1) EMR-1 is not 
only restricted to the NE but can diffuse into the nuclear 
interior; (2) EMR-1 localizes to patches of the NE devoid 
of LMN-1; or (3) EMR-1 is recruiting proteins to these 
DNA elements that inhibit access of LMN-1. 

We consider the first possibility unlikely because 
EMR-1, a transmembrane domain-containing protein. 
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Figure 3 EMR-1 associates with muscle and neuronal genes. (A) Comparison of the regions enriclied for LMN-1 and EMR-1 in cliromosomes I, III 
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'EMR-1 only' elements as defined in Materials and methods. (C) Genome browser views of LMN-1 and EMR-1 occupancy in regions containing 'EMR-1 
only' elements (grey boxes) at promoters (left panel) and in intergenic regions (right panel). (D) Distribution of 'EMR-1 only' probes in promoters, introns 
and exons relative to the genome. Probability values from one-sided binomial tests are indicated as in (B). (E) Gene Ontology (GO) analysis of genes 
included in 'EMR only' elements. GO terms with P< 0.05 were considered significant. Numbers in the pie chart indicate the number of genes in each 
category. (F) LMN-1 and EMR-1 occupancy for various gene sets defined by expression data. Number of genes: intestine (555), muscle (359), serpentine 
(325), neuron (14), germ line (618) and ubiquitous (174). Serpentine genes encode chemoreceptors expressed in a few neurons only. Boxplots (as in 
Figure 1) show the range of MA2C scores averaged over gene bodies. Probability values from Wilcoxon rank sum tests are indicated: ""PkO.OS, """"Pk 10"^. 
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had never been visualized in the nuclear interior [37,38]. 
Although NE accumulation of EMR-1 is LMN-1- 
dependent [36], we cannot rule out the second possibil- 
ity since biochemical purification of emerin from human 
HeLa cells identified protein complexes that lack lamins 
[39]. Moreover, small LMN-l-free patches of the NE 
involved in anchoring of 'EMR-1 only' elements would 
typically be undetectable at the resolution of light mi- 
croscopy. To clarify whether EMR-1 could recruit add- 
itional factors to the 'EMR-1 only elements, thereby 
sterically inhibiting Dam::LMN-l access, we expressed 
Dam::LMN-l in mutants homozygous for the emr-1 null 
allele ^/d 79 [30]. If this hypothesis is true, in the absence 
of EMR-1, these elements should be able to interact with 
LMN-1. Based on the redundancy between EMR-1 and 
LEM-2 during embryogenesis [29] we also included the 
lem-2{tml582) allele in our studies. Consistent with the 
absence of obvious growth phenotypes, mutations of 
emr-1 or did not produce large-scale changes in 

LMN-1 distribution (Pearson correlation of 0.94 and 
0.89 between wild type and emr-l and lem'2 mutants, 
respectively; Figure 4A; Figure S7 in Additional file 1). 
However, in nematodes lacking EMR-1 or LEM-2, 
the ability of LMN-1 to interact with chromatin was 
mildly but significantly reduced, as can be seen by the 
lower MA2C score observed for LADs in these mutants 
(Figure 4B). In contrast, in both mutants, LMN-1 was 
localized more strongly to areas that were designated as 
'EMR-1 only in wild- type nematodes (Figure 4B, C). 
This effect was specific for 'EMR-1 only' elements since 
'LMN-1 only' elements behaved like the rest of the 
LADs, with a mild and general decrease in association 
with the NE (Figure 4B). These data support the idea 
that EMR-1 (and LEM-2) could recruit transcription or 
chromatin remodeling factors that inhibit LMN-1 asso- 
ciation with a subset of 'EMR-1 only' genomic locations. 

Although lack of EMR-1 expression did not produce 
massive reorganization of the genome, we identified re- 
gions that were detached from the nuclear periphery upon 
EMR-1 deletion (Figure 4C; Table SI in Additional file 2). 
These regions were specifically enriched in introns 
(Figure 4D) and GO analyses revealed that they were asso- 
ciated with genes involved in development of reproductive 
organs and muscles (Table S6 in Additional file 3). 

Depletion of EMR-1 and LEM-2 increases gene expression 
globally but reduces transcription of chromatin modifying 
enzymes 

The results above suggest that EMR-1 and LEM-2 may 
regulate the expression of certain gene classes. To test 
this directly, we performed RNA sequencing in wild- 
type animals and in animals depleted for EMR-1 and/or 
LEM-2. Data from wild-type nematodes confirmed that 
EMR-1 was associated with silent and poorly transcribed 



regions (Figure 5A). Similar results were obtained for 
LADs (Figure S8A in Additional file 1), supporting the 
notion of the nuclear periphery being a zone of generally 
repressed chromatin. Genes included in 'EMR-1 only' el- 
ements also showed reduced expression relative to genes 
in gaps. However, expression within 'EMR-1 only' re- 
gions was slightly but significantly higher than the ex- 
pression across all EADs (Figure 5A). 

Since LMN-1 DamID profiles of emr-1 or lem-2 single 
mutants did not reveal dramatic changes in LMN-1 do- 
main organization, we analyzed gene expression in ani- 
mals lacking expression of both EMR-1 and LEM-2. To 
bypass the embryonic lethality produced by the simul- 
taneous deletion of EMR-1 and LEM-2, we knocked 
down emr-1 expression by RNA interference (RNAi) in 
lem-2 mutants. LI larvae homozygous for the lem-2 
(tmlS82) allele were fed with bacteria expressing emr-1 
double-stranded RNA. These nematodes developed to 
adulthood at the normal rate but produced a large num- 
ber of dead embryos, which verified the effectiveness of 
the RNAi treatment [29]. As illustrated in Figure 5B, a 
general increase in RNA levels for all genes was ob- 
served in both singly and doubly depleted animals: in 
lem-2(tml582); emr-l(RNAi) doubly depleted animals, 
the RNA levels of more genes were significantly altered 
relative to wild type (9.7% of protein-coding genes com- 
pared with 2.2% and 0.48% in emr-l(gkll9) and lem-2 
(tml582) single mutants, respectively; Figure 5B; Figure 
SBC in Additional file 1), most likely reflecting the re- 
dundant functions of EMR-1 and LEM-2. 

Analysis of the genes that lose their NE interaction in 
emr-1 mutants revealed a negative correlation between 
NE association and gene expression (Figure S8E in 
Additional file 1). Combined with the decreased associ- 
ation of chromatin to the nuclear periphery observed 
by LMN-1 DamID in the single mutants (Figure 4B), 
this result clearly indicates a role for EMR-1 and LEM-2 
in gene repression. However, consistent with previous 
studies [34,40], we observed that not all the genes that 
were up-regulated in the emr-1 mutant changed their 
position relative to the NE (Figure S8D in Additional 
file 1). This indicates that multiple processes influence 
the relationship between NE association and gene ex- 
pression (see Discussion). 

When we analyzed genes expressed in specific tissues, 
we also found that effects were maximal in the doubly 
depleted animals. Interestingly, an increase in expression 
of intestinal, muscle and serpentine genes was observed 
whereas germ line and ubiquitous genes were repressed 
(Figure 5C). Because the number of neuronal genes in 
this data set was very low (n = 14) we expanded our ana- 
lysis to include 338 genes identified as pan-neuronal 
genes in larvae [41] (Table S2 in Additional file 2). 
Of these, 141 genes are associated with the nuclear 



Gonzalez-Aguilera et al. Genome Biology 2014, 15:R21 
http://genonnebiology.conn/201 4/1 5/2/R21 



Page 10 of 19 



WT emM (gk119) lem-2 (tm1582) 



1.4. 

CD 
O 

o 
CO 

O 0. 

C\J 

< 



1.2 










Chrlll 







ChrlV 



B 



LMN-1 Domains 



o 
o 

0) 

O 

C\J 

< 



2 

1.5 

1 ' 
0.5^ 



I I * 



I- 

5)^ 



Dam::LMN-1 



EMR-1 Only 



LMN-1 Only 




LMN-1 
EMR-1 



EMR-1 Only 



LMN-1 2.5 
emr-1A ^ 

2.5 



-2 
2.5 

-2 



4kb 



4.5kb 



LMN-1 2.5- 
emr-1A _2_ 
2.5 



LMN-1 



EMR-1 



-2 
2.5 

-2 



Detached regions 





8.5kb 



4.7kb 



D 



^ 60- 
oS 50- 
^ 40- 
f 30- 
^ 20 
10 



I Genome 

I Detached regions 



^^ri *** I 

M 

Promoters exons introns 
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Figure 5 Mutations in emr-1 and lem-2 de-repress gene expression. (A) Boxplot (as in Figure 1) of expression levels of genes included in 
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of all genes in wild type {WT), Iem-2(tml582), emr-l(gkl 19) and emr-l(RNAi) Iem-2(tml582) ('double') animals. (C) Expression of tissue-specific genes 
grouped as described in Figure 3. ""Pk 10"^ """"P < 10"^ and ***p< 10"". (D) Gene Ontology analysis of genes up- (green) or down-regulated (black) 
in lem-2 (tm 1582), emr-l(gkl 19) and emr-l(RNAi) lem-2(tml 582) animals compared to wild type (P<0.01). (E) Chromatin modif/ing enzymes with 
altered expression in emr-1 (RNAi) Iem-2(tml582) nematodes. Most genes are down-regulated, except underlined genes, which are up-regulated. 
Genes also affected in the emr-1 {gkl 19) single mutant are shown in bold. 



Gonzalez-Aguilera et at. Genome Biology 2014, 15:R21 
http://genonnebiology.conn/201 4/1 5/2/R21 



Page 1 2 of 1 9 



periphery as judged from positive EMR-1 DamID MA2C 
scores. Interestingly, the pan-neuronal genes were signifi- 
cantly up-regulated in doubly depleted animals (Figure 
S6C in Additional file 1). These results were confirmed 
by GO analyses since GO categories related to the ner- 
vous system and locomotion were enriched among the 
significantly up-regulated genes (Figure 5D; Table S8 in 
Additional file 3). In contrast, down-regulated genes 
were associated with GO terms related to reproduction 
and basic cellular processes like cell cycle or transcription 
(Figure 5D). Intriguingly, genes related to chromatin 
organization were over-represented among the transcripts 
with reduced expression {P = 1.17e-6). In fact, the expres- 
sion of 18 of the 38 chromatin modifying enzymes de- 
scribed in C elegans [42] was significantly affected in the 
doubly depleted animals (47%, compared to 9.7% of all 
genes; Figure 5E; Table S8 in Additional file 3). We cannot 
rule out the possibility that some of the changes in expres- 
sion levels of chromatin modifier enzymes could be sec- 
ondary effects of chromatin reorganization upon EMR-1 
and LEM-2 depletion. Nevertheless, these data support 
the idea that EMR-1 and LEM-2 could be especially im- 
portant for regulating the activity of genes encoding 
chromatin-remodeling activities. 

EMR-1 Is required for proper neuromuscular junction activity 

Our results are in concordance with a role of emerin/ 
EMR-1 and LEMD2/LEM-2 in muscle formation and 
activity as has been reported in various organisms 
[30,31,43]. However, to our knowledge, neither emerin 
nor LEMD2 has been implicated in neuronal processes. 
Since both NE tethering and expression of genes re- 
lated to the nervous system and synaptic transmission 
were altered in emr-1 and lem-2 mutants (Figures 3, 4 
and 5), we explored this possibility. Indeed, among the 
genes related to synaptic transmission, the expression of 
those encoding neurotransmitter vesicle factors, chan- 
nels, and acetylcholine receptor proteins was affected 
(Figure 6A). To investigate in vivo whether EMR-1 and 
LEM-2 play a role in neural physiology, we performed 
aldicarb sensitivity experiments to assess synaptic trans- 
mission activity. Aldicarb is a chemical analogue of the 
neurotransmitter acetylcholine and inhibits the activity of 
cholinesterase at neuromuscular synapses. Consequently, 
acetylcholine accumulates in the synaptic cleft, causing 
continuous activation of muscular acetylcholine receptors, 
muscle hyper-contraction and paralysis [44]. Alteration in 
aldicarb sensitivity can therefore be used to detect changes 
in acetylcholine release and perception. 

We determined the time-course for aldicarb-induced 
paralysis in wild type, emr-1 and lem-2 single mutants 
and in lem-2 nematodes heterozygous for the emr-1 mu- 
tation {emr-l(gkll9)h; lem-2(tml582)). We did not test 
for EMR-1 /LEM-2 double depletion because RNAi is 




Time (min) 

Figure 6 Synaptic function is altered in emr-1 mutants. (A) Fold 
increase in expression of selected genes involved in synaptic 
transmission in emr-1 (RNAi) Iem-2(tml582) animals. Error bars indicate 
minimum and maximum values. (B) Aldicarb sensitivity experiments 
in wild type [WT), emr-l(gkl 19), Iem-2(tml582) and emr-1 (gkl 19)/+ 
Iem-2(tml582) mutants. Synaptic function was evaluated based on 
paralysis induced by the cholinesterase inhibitor aldicarb. Data 
points show the mean of four experiments and error bars indicate 
standard error of the mean. 

V J 

less effective in C. elegans neurons and the homozygous 
double mutants proceeding from heterozygous animals 
arrest during larval development [30]. Compared with 
wild-type animals, emr-1 mutants were hypersensitive to 
aldicarb (50% of animals were paralyzed after approxi- 
mately 35 minutes for emr-1 mutants versus approxi- 
mately 55 minutes for wild- type animals), whereas lem-2 
mutants retain normal aldicarb sensitivity (Figure 6B). 
Heterozygous emr-l(gkll9)/-\- mutants were equally sen- 
sitive to aldicarb as wild-type animals, indicating that a 
single copy of emr-1 is sufficient to maintain proper 
neuromuscular junction activity. 

At this moment, we cannot distinguish whether the 
paralysis observed in emr-1 mutants is due to changes in 
neurotransmitter release or muscle perception, or a 
combination of the two. Interestingly, lem-2 mutants 
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present alterations in muscle structure and function [30] 
and, unlike emr-l mutants, do not show altered aldicarb 
sensitivity. This indicates that EMR-1 and LEM-2 are 
not redundant in all their functions, making C. elegans a 
good model to study the tissue-specific roles of these 
two proteins. 

Discussion 

A comparison of the genome-wide DNA association 
profiles of LMN-1, the only lamin protein in C. elegans, 
and of the INM protein emerin/EMR-1 indicates that 
EMR-1 defines two types of chromatin domains: a major 
type consisting of LADs and a minor type devoid of 
LMN-1 and enriched for genes expressed in specific 
tissues. LADs are conserved from embryos to adults 
and epigenetic marks suggest they are transcriptionally 
silent at both developmental stages (this study and [34]). 
In contrast, genes associated with EMR-1 only show 
higher expression and belong to functional categories 
that respond to stimuli. The difference between these 
two classes of chromatin domains might be mediated by 
posttranslational modifications of EMR-1 as well as 
other NE associated proteins. For instance, phosphoryl- 
ation of human emerin has been shown to be important 
for the association with BAP [45,46]. 

The discovery of regions associated with EMR-1 but 
not with LMN-1 was surprising, since the two proteins 
physically interact at the NE and LMN-1 is required to 
localize EMR-1 properly [36]. Interestingly, although we 
could find 'EMR-1 only' elements in all chromosomes, 
they were especially abundant in the X chromosome. 
We do not know why the X chromosome has fewer in- 
teractions with LMN-1 compared to the autosomes, nor 
can we explain the particular enrichment of 'EMR-1 
only elements on the X chromosome. Our experiments 
were performed with C. elegans hermaphrodites, in 
which both X chromosomes undergoes approximately 
two-fold chromosome-wide transcriptional repression to 
achieve dosage compensation [47]. It is possible that the 
interaction with EMR-1 or other INM proteins could be 
important for this transcriptional repression. Our data 
suggest that the anchoring of the X chromosome to the 
NE could be more pronounced than previously inferred 
from LEM-2 ChIP data [18] or by LMN-1 DamID ([34] 
and this study). 

We found that 'EMR-1 only elements were sur- 
rounded by chromatin interacting with both LMN-1 and 
EMR-1. This indicates that 'EMR-1 only' elements are 
located at or very close to the NE, but raises the ques- 
tion as to why LMN-1 is not able to bind to these chro- 
matin regions. Given that the absence of EMR-1 and/or 
LEM-2 allows the 'EMR-1 only' elements to interact 
more frequently with LMN-1, we hypothesize that EMR- 
1 and LEM-2 physically prevent LMN-1 from binding to 



these regions. It is known that human and mouse 
emerin interacts physically with transcription factors in- 
volved in muscle development, including p-catenin [23] 
and Lmo7 [24], and with chromatin associated factors, 
such as BAP and HDAC3 [13,48]. Similarly, the mamma- 
lian LEM-2-like protein MANl interacts with Smads at 
the NE to regulate transforming growth factor- 13 signal- 
ing [49]. C. elegans EMR-1 and LEM-2 may also associ- 
ate with DNA-binding proteins and such protein 
complexes may hinder LMN-1 from binding to 'EMR-1 
only' elements. 

The NE is generally associated with repressed DNA 
and several of our observations suggest that EMR-1 and 
LEM-2 may also be involved in anchoring of silent het- 
erochromatin to the nuclear periphery in specific tissues. 
Pirst, genes within EADs are less expressed than genes 
in gaps. Second, we have demonstrated that, during 
adulthood, depletion of EMR-1 and/or LEM-2 caused a 
partial release of LADs from the NE and increased tran- 
scription of several hundred genes. These results are 
consistent with the delocalization and overexpression of 
repetitive heterochromatic transgene arrays observed 
upon simultaneous depletion of EMR-1 and LEM-2 in 
adult nematodes [50]. Interestingly, although chromatin 
associated with EMR-1 is less transcribed than chroma- 
tin found in gaps, genes in 'EMR-1 only' elements are 
significantly more expressed than genes in other EMR-1 
domains. It is possible that the accumulation of tran- 
scription factors around 'EMR-1 only' elements might 
poise these regions for a rapid response to environmen- 
tal stimuli. Supporting this idea, 'EMR-1 only' elements 
are enriched for GO categories related to signaling path- 
ways and transcription. However, we have also noticed 
that gain or loss of NE anchoring do not always correl- 
ate with changes in gene expression. This has also been 
observed in other organisms. Por instance, approxi- 
mately one-third of genes that move towards or away 
from the NE during mouse embryonic stem cell (ESC) 
differentiation do not change transcriptional activity 
[40]. Similarly, changes in expression of genes that gain 
or lose lamin interaction after ecdysone treatment is 
marginal in Drosophila Kc cells [19]. Pinally, mouse 
ESCs lacking B type lamins show altered transcription 
profiles but the alterations do not correlate with lamin 
Bl association in wild-type ESCs [51]. The reason for 
this lack of concordance between NE anchoring and si- 
lencing is not clear. Peric-Hupkes and colleagues [40] 
proposed that genes that move away from the NE dur- 
ing ESC differentiation are unlocked' to become acti- 
vated at later steps of development. Likewise, results 
from a recent study in human adipose stem cells, indi- 
cate that dissociation of promoters from the nuclear 
lamins produces chromatin changes that enhance tran- 
scriptional permissiveness but do not necessarily elicit 
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gene activation [52]. The molecular details of these 
intermediate states are unknown. The observation in 
C. elegans that the sequential action of two H3K9 
methyltransferases promotes first anchoring of DNA to 
the NE and afterwards transcriptional repression has 
provided mechanistic insight into the uncoupling of the 
two processes [34]. 

Our data also indicate that genes specifically expressed 
in muscles and neurons are enriched in EADs. Further- 
more, expression levels of genes related to locomotion 
and behavior were affected in emr-1 and lem-2 mutants. 
Several studies on Emery-Dreifuss muscular dystrophy 
patients and model organisms support the role of emerin 
and LEMD2 in muscle integrity, but their potential role 
in neurons had not been explored. Here we show that in 
C. elegans, localization and expression of genes related 
to neurotransmitter recycling, acetylcholine receptor 
network and ion-channels are affected in emr-1 and 
lem-2 mutants. Interestingly, neuromuscular junction 
activity is strongly affected in emr-1 but not in lem-2 
mutants, indicating that the two genes are not redun- 
dant in all their functions. The role of emerin in neur- 
onal activity may not be restricted to C. elegans since 
emerin knockout mice suffer relatively mild muscle 
defects and altered motor coordination [8]. In rodents, 
expression of emerin but not LEMD2 has been detected 
in neuroretinal cells, while both proteins are equally 
present in muscles [53]. Moreover, mutations in several 
other NE components, such as lamin A, BAP, lamin Bl, 
and the LINC complex (which connects the NE to the 
cytoskeleton), are associated with neuropathies affecting 
both the peripheral and central nervous system [54-58]. 
Future analyses are needed to investigate whether 
emerin also contributes to proper neuronal physiology 
in humans. 

Finally, we have demonstrated that, in adult nema- 
todes, the NE-associated heterochromatin domains are 
highly enriched at the ends of autosome arms and at 
the left arm of chromosome X. This configuration is 
very similar to the one observed during embryogenesis 
[18,34]. However, we also found differences in NE- 
genome associations during development. These include 
few local changes in LADs placed in the chromosome 
arms and a general decrease in association of LMN-1 
with the central part of the chromosomes in adults. 
Developmental changes in chromatin distribution have 
been observed in several organisms. The location of in- 
tegrated small transgene arrays has been reported to be 
developmentally regulated, being randomly positioned in 
embryonic nuclei and accumulated in the nuclear inter- 
ior in larvae [59]. The integration sites of these trans- 
gene arrays are unknown but their repositioning away 
from the NE coincides with transcriptional activation 
[59]. In murine embryonic stem cells, approximately 



10% of genes change their interaction with the nuclear 
lamina during in vitro differentiation [40]. Together, 
these results suggest that, across evolution, global nu- 
clear organization is established early in development 
and only a few specific regions change their localization 
in response to developmental requirements or changes. 

Conclusions 

The comparison of chromatin associated with the con- 
served NE proteins LMN-1 /lamin and EMR-1 /emerin 
presented here reveals a large overlap between LMN-1- 
and EMR-1 -associated domains but also the existence of 
numerous elements that interact with EMR-1 only. 
Muscle and neuronal genes are enriched in these 'EMR- 
1 only elements and are deregulated in animals depleted 
for EMR-1 and LEM-2. Moreover, differential expression 
of EMR-1 and LEM-2 was observed across tissues. 

The requirement for EMR-1 for the proper spatial and 
transcriptional control of muscle and neuronal genes 
demonstrated here prompts the idea that variation in 
NE composition among tissues may contribute to the 
correct adaptation of the organism to the environment. 
Combined with the defects in neuromuscular junction 
activity in emr-1 mutants reported here, we propose that 
development of Emery-Dreifuss muscular dystrophy may 
involve altered nuclear architecture and gene expression, 
especially in muscles and neurons. 

Materials and methods 

Strains and plasmids 

Strains carrying a single-copy insertion of Dam-fusion 
constructs were made using the MosSCI technique 
[60]. Strain and plasmid construction are described in 
Additional file 4. 

Immunofluorescence and western blot 

Nematodes grown at 20°C were heat-shocked 1 h at 33°C 
and left to recover for 2 h at 20°C. Embryos were obtained 
and processed as described [61]. Anti-MYC primary anti- 
body (Sigma-Aldrich, catalog number C3956, St. Louis, 
MO, USA) was diluted 1:100 for immunofluorescence and 
1:400 for western blotting. Secondary goat anti-rabbit 
Alexa Fluor 633-conjugated antibody (Invitrogen, catalog 
number A21071, Carlsbad, CA, USA) diluted 1:1,000 and 
peroxidase-conjugated antibody (Sigma-Aldrich, catalog 
number A0545) diluted 1:5,000 were used for immuno- 
fluorescence and western blotting, respectively. 

DamID experiments 

DamID experiments were performed using three bio- 
logical replicates and all cultures were grown in parallel. 
DamID strains were synchronized from eggs prepared 
by hypochlorite treatment of gravid adults fed with 
OP50 E, coli and left to hatch overnight at 20°C in 
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S -medium. Then, from each strain, approximately 
35,000 Lis were grown in 50 ml S -medium containing 
GM119, a Dam" E, coli strain (F' dam-3 dcm-6 metBl 
galk2 galT22 lacYl tonASl tsx-78 supE44 mtl-1). This 
food source was used to avoid contamination by methyl- 
ated E, coli DNA. Cultures were grown with continuous 
agitation (180 rpm) at 20°C for 53 h. At this time 
point the cultures were highly enriched for non-gravid 
young adults. Nematodes were frozen at -80°C until 
further processing. 

Methylated genomic DNA (gDNA) was purified and 
amplified from 30 mg nematodes using a protocol based 
on [62] with some modifications. Briefly, 2.5 (ig of gDNA 
isolated with DNAeasy kit (QIAGEN, Venlo, Limburg, 
Netherlands), was digested with 10 units Dpnl (New 
England Biolabs, Ipswich, MA, USA) overnight in 10 \A 
to cut at G^'ATC sites. Dpnl was inactivated (80°C, 20 
minutes) and then gDNA was ligated to double-stranded 
adaptors with 5 units T4 DNA ligase (5 U/(il; Roche, 
Basel, Switzerland) in a volume of 20 (il. After inactiva- 
tion of the ligase (65°C, 10 minutes), DNA fragments 
were digested with 5 units Dpnll (New England Biolabs) 
in a final volume of 80 \A to cut non-methylated GATC 
sites, thereby preventing PGR amplification of non- 
methylated gDNA. Methylated DNA was amplified using 
adaptor-specific primers. PGR reactions were performed 
using 1 \A PGR Advantage enzyme mix (50x; Glontech, 
Otsu, Shiga, Japan), 10 \A of Dpnll digestion and the fol- 
lowing program: (1) 68°G for 10 minutes, (2) 94°G for 1 
minute, (3) 65°G for 5 minutes, (4) 68°G for 15 minutes, 
(5) 94°G for 1 minute, (6) 65°G for 1 minute, (7) 68°G 
for 10 minutes, (8) go to step 5 3x, (9) 94°G for 1 mi- 
nute, (10) 65°G for 1 minute, (11) 68°G for 2 minutes, 
(12) go to step 9 20 X. To achieve sufficient material, two 
PGR reactions from the same Dpnll digestion were 
pooled and purified with Qiaquick PGR purification kit 
(QIAGEN). Amplified Dam-methylated DNA was la- 
beled and hybridized by the Roche Nimblegen Service 
Laboratory. Dye orientation of experiments is indicated 
in Table S4 in Additional file 2. 

DamID data processing 

Nimblegen 2.1 M whole genome tiling arrays, with 50 
bp probes, designed against WS180 (ce5) genome as- 
sembly, were used for all experiments. Data were ana- 
lyzed with MA2G [33], which was used to normalize the 
log2 ratio (log2 [Dam::fusion/GFP::Dam control]) of each 
probe based on the probe behavior estimated by its GG 
content, and then to smooth the value by assigning 
the median across sliding windows of 300 bp. The re- 
sultant values are MA2G scores. Subsequently, data were 
subjected to quantile normalization over all available 
replicates and backgrounds to facilitate comparison be- 
tween strains. Finally, the mean of normalized MA2G 



scores from the three independent biological replicates 
was calculated. We then converted the chromosome co- 
ordinates to WS190 (ce6) genome assembly using the 
Lift-Over tool from the UGSG genome browser [63] and 
used these values for the subsequent domain calling and 
sliding window analyses. In some cases, to facilitate 
chromosome-scale data visualization by reducing the 
number of data points, MA2G scores were averaged 
within non-overlapping 100 kb windows across the gen- 
ome as indicated in the figure legends. 

Sliding window analysis was performed using the left 
boundaries corresponding to the left edges of all EMR-1 
domains (Table SI in Additional file 2). Both sides of 
boundaries and sliding window analyses with LMN-1 
boundaries provided the same results. Boxplots for gene 
classes were generated by averaging MA2G scores of 
probes included in regions between the transcription 
start site and end site. Goordinates for these sites were 
obtained from the UGSG genome browser (WS190 gen- 
ome assembly). Analyses for enrichment at promoters 
(up to 3 kb before the transcription start site), exons and 
introns were performed using the GEAS program [64]. 
GhlP-chip data of young adult H3K4me3 and H3K36me3 
marks [35], and the embryonic LEM-2 mark [18] are from 
modENGODE [65] (accession IDs modENGODE_3552, 
modENGODE_3559, modENGODE_2729, respectively). 
Embryonic LMN-1 DamID data [34] were obtained 
from the Gene Expression Omnibus (GEO; accession 
GSE37226). Gene sets used in tissue-specific analyses were 
defined previously and are listed in Table S2 in Additional 
file 2: muscle [66], intestine [67], ubiquitous, serpentine, 
germ line [68], and neuron [69]. 

Definition of LMN-1- and EMR-1 -associated domains 

LMN-1- and EMR-1 -associated domains were defined 
as previously described [18]. Briefly, we transformed 
positive and negative MA2G scores to +1 and -1, re- 
spectively. We averaged the binary values in 200-probe 
windows (approximately 10 kb), sliding one probe (50 
bp offset) across the genome to identify windows with 
high average binary values. To have an estimation of the 
false positive windows, we performed the same proced- 
ure on a GFP::Dam/GFP::Dam MA2G normalized data- 
set. From this, we defined positive windows when a 
binary window value is >0.8, at which the false discovery 
rate is <5% (Figure S2 in Additional file 1). We then 
joined any overlapping (>1 bp) windows to generate the 
LMN-1- and EMR-1 -associated domains (Table SI in 
Additional file 2). 

The definition of chromosomal arms and chromo- 
somal centers for comparison of embryonic and adult 
LMN-1 association profiles (Figure 2D) was based on 
the boundaries of LADs in adult nematodes. Goordinates 
of chromosome arms: chromosome I, 1 to 3,745,632 and 
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10,809,938 to 15,072,421; chromosome II, 1 to 4,708,341 
and 11,877,168 to 15,279,323; chromosome III, 1 to 
3,508,994 and 9,947,268 to 13,783,681; chromosome IV, 
1 to 7,317,812 and 12176625 to 17,493,785; chromosome 
X, 1 to 4,191,936. 

To define 'EMR-1 only' elements, we selected probes 
with a MA2C score >1 in EMR-1 DamID experiments 
and <0 in LMN-1 DamID experiments. We defined 
'EMR-1 only elements as the regions that contain at 
least 10 positive probes, each separated <500 bp from 
the neighboring positive probe (Table SI in Additional 
file 2). We performed the same procedure to identify 
'LMN-1 only elements (Table SI in Additional file 2). 

RNAi and RNA extraction 

To simultaneously deplete EMR-1 and LEM-2, BN19 ani- 
mals homozygous for the lem-2(tml582) allele were 
subjected to emr-l RNAi-mediated knockdown. RNAi 
feeding bacteria were grown in LB at 37°C with agitation 
(180 rpm) until ODeoo = 0-8 and then induced with 1 mM 
isopropyl |3-D-l-thiogalactopyranoside (IPTG) during 19 
h at 20°C. The RNAi bacteria expressed double-stranded 
RNA corresponding to either full-length emr-l cDNA or, 
as negative control, 185 bp of unrelated sequence from 
the L4440 empty plasmid. The next day, the bacteria 
were collected and resuspended in 1/20 initial volume 
of S-medium with 1 mM IPTG. Approximately 40,000 
synchronized LI nematodes obtained by hypochlorite 
treatment and hatching overnight at 20°C without food 
were grown in 50 ml of S-medium containing 5 ml of 
resuspended RNAi bacteria, 100 (ig/ml ampicillin and 1 
mM IPTG at 20°C during 53 h with continuous agita- 
tion (180 rpm). At this point the cultures accumulated 
non-gravid young adults. Nematodes were recovered, 
frozen immediately in liquid N2 and kept at -80°C until 
further processing. Part of the culture was grown for an 
additional day, which confirmed the expected embry- 
onic lethal phenotype in the EMR-1 LEM-2 double de- 
pletion. RNA extraction was performed using 20 to 30 
mg of nematodes. The nematode cuticle was ruptured 
by three rounds of freezing in liquid N2 and thawing at 
30°C and RNA was purified using the RNeasy Mini kit 
(QIAGEN). 

Expression profiling 

Sequencing libraries were prepared using the TruSeq 
RNA sample preparation kit according to the manufac- 
turers instructions (Illumina, San Diego, CA, USA). 
Briefly, polyA-containing mRNA was isolated from 1 [ig 
of total RNA using poly-T oligo- attached magnetic beads 
and processed as described [70]. Pools of three indexed 
libraries were mixed (multiplexed) at equimolar ratios 
to yield a total oligonucleotide mix of 10 nM. Finally, 
the resulting libraries were sequenced on the Illumina 



Genome Analyzer IIx platform to generate 150 bp 
single-end reads. Three pooled indexed libraries were se- 
quenced in each flow cell lane. 

Raw RNA sequence data were processed [70] and 
aligned with the WS231 genome version using TopHat. 
Next, results were sorted and indexed using Samtools. 
Coverage vectors for each sample were extracted under 
R statistical environment using Rsamtools library. The 
C. elegans GTF annotation from Ensembl was applied as 
reference for the exon genomic ranges, parsed into R data 
structures. Exon expression was calculated as the median 
coverage for every exon genomic range as stated by anno- 
tation. Exon expression data from our nine different 
samples (three N2, two BN19, two BN20, and two BN24) 
were normalized by quantile normalization and locus 
expression was defined as the median of all the expression 
hits from distinct exons annotated to a single locus. 
Library siggenes in R was applied to select the differen- 
tially expressed loci using the false discovery rate-based 
SAM method [71]. Genes presenting a A value >4 were 
considered to be significantly altered (Table S3 in 
Additional file 2). Resulting data are available with ac- 
cession numbers summarized in Table S4 in Additional 
file 2. For the analysis of emr-l (gkl 19) and lem-2 
(tml582) mutants, the sets of genes expressed in differ- 
ent tissues were the same as in the DamID analysis 
(Table S2 in Additional file 2). 

Gene Ontology analysis 

GO analyses were performed using DAVID [72] and only 
the 'Biological Processes' tree was used in our study. For 
the study of DamID 'EMR-1 only' elements, as they were 
enriched in promoters and regulatory regions, the clos- 
est but not overlapping genes were selected. In the other 
cases, we selected genes with >50% of the open reading 
frame included within the domains. 

Aldicarb sensitivity experiments 

Sensitivity to aldicarb was determined by assaying the 
time course of the onset of paralysis following acute ex- 
posure of a population of animals to this drug. In each 
experiment 20 to 25 young adults were placed on freshly 
made drug plates (1 mM aldicarb; Fluka Sigma- Aldrich 
ref. 33386) and prodded every 10 minutes over a 2 h 
period to determine if they retained the ability to move. 
Nematodes that failed to respond to touch were classi- 
fied as paralyzed. Each experiment was repeated four 
times by two independent experimenters. 

Data access 

DamID and RNA sequencing data performed in this study 
are available from GEO [73] accession IDs GSE44188 
and GSE44682, respectively (described in Table S4 in 
Additional file 2). 
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Additional file 1: Figure SI. DamID fusion proteins localize properly at 
the NE. Figure 52: quality assessment of LMN-1 and EMR-1 DamID data. 
Figure S3: LMN-1 and EMR-1 associate with the ends of chromosomes 
and silenced DNA. Figure S4: analysis of 'EMR-1 only' and 'LMN-1 only' 
elements. Figure S5: differential expression of EMR-1 and LEM-2. Figure 
S6: genes expressed in different tissues show distinct association with 
EMR-1 and LMN-1. Figure S7: analysis of LMN-1 DNA association in wild 
type and emr-1 and lem-2 mutants. Figure S8: expression analyses of 
genes associated with LMN-1 and EMR-1. 

Additional file 2: Table SI. The coordinates for EADs, LADs, 'EMR-1 

only' elements and 'LMN-1 only' elements. Table S2: genes that are 
expressed in specific tissues. Table S3: RNA sequencing data. Table S4: 
GEO accession numbers. 

Additional file 3: Tables S5 to S10: Information on significant GO 
categories. 

Additional file 4: Supplementary Materials and methods. 
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